Expansion in the "distance" from Hc2(T) line for the mixed state of BCS 

superconductors 
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We develop a description of the mixed state of type II superconductivity vaUd within a wide range 
of temperatures and external magnetic fields. It is based on the quasiclassical version of microscopic 
BCS theory and employs an expansion in the "distance" from the Hc2{T) line. We prove that in 
a clean metal with spherical Fermi surface and isotropic pairing interaction the superconducting 
condensation always produces a hexagonal vortex lattice. 



The theoretical description of the mixed state is a dif- 
ficult problem: the equations are nonlinear, the solution 
for the ground state is inhomogeneous and selfconsistency 
is required. Analytical progress is facilitated if a small 
parameter is present in the theory. In particular, a per- 
turbative treatment partially removes the need for self- 
consistency. Near the second order phase transition the 
order parameter is small. The Ginzburg-Landau (GL) 
approach is a tool to describe such situations, and it was 
proven successful on numerous occasions, including the 
very discovery of the upper critical field Hc2 and the vor- 
tex nature of the mixed statil. The Hc2{T) line can 
be calculated from the microscopic BCS theory|-.Qf su- 
perconductivity: in some sjopmetric cases exactlyo'Q and 
more often approximatelyQ'0. The framework here is the 
Gorkov formalism or, even more conveniently, the quasi- 
classical Eilenberger equations. 

The aim of this paper is to develop a rigorous pertur- 
bative method to describe the mixed state. An analytical 
description would allow to follow general trends as one 
moves onT — H plane. To find the structure of a vortex 
solid at an arbitrary (T, H) point directly is only possi. 
by numerical means and constitute demanding wori 
Therefore, one often starts from the Hc2{T) line and ex- 
ploit in some way the smallnesa.of the pairing amplitude 
in its vicinity (see for examplea). It is worth emphasiz- 
ing, however, that the order parameter is a function of 
spatial coordinates while to perform an expansion un- 
ambiguously a small number is required. Below we argue 
that the "distance" to the Hc2 [T) line can be used as such 
a small number and we demonstrate how to develop the 
expansion. We find that this parameter remains small 
on a very large portion of the superconducting region of 
T — H plane, indicating at possible large range of validity 
of the expansion. The idea to use the "distance" to the. 
Hc2{T) line was recently introduced by B. RosensteinEI 
for the analysis of fluctuations in the framework of GL 
theory. For the sake of simplicity we consider in this pa- 
per the case of a clean metal with spherical Fermi surface 
(FS) and isotropic pairing interaction. As an immediate 
result we are able to prove that superconducting con- 
densation produces the hexagonal vortex lattice at any 
temperature. A further development of this formalism to 
more complicated cases is possible. It could provide, in 
particular, valuable information on vortex lattices in un- 



conventional superconductors^ which are under scrutiny 
at present, but have been analyzed mostly in the GL 
framework so far. 

The Eilenberger equations relate the quasiclassical 
Green functions, both normal (7(a;n;v,x) and anomalous 
/(a;„;v,x), /(cj™ju_x), and tiie superconducting order 
parameter A(x)tilll3 (see Ref.ll3 for review). Using the 
notations v for the unit vector in the direction of the 
Fermi velocity, a;„ = 2TTt{n +1/2) for Matsubara fre- 
quencies and Di = di — ie*Ai for covariant derivatives, 
they can be written as 



(w„ + vD/2)/ = Ag, 
(w„-vD72)/ = A*g, 



(1) 

(2) 
(3) 



The last relation is the constraint on the quasiclassical 
Green functions which can be used instead of the differen- 
tial equation for g. Below we use the following units: Tc 
for temperature {t — T/Tc is reduced temperature) and 
the gap function, R = Tivp/Tc for distances and Hn = 
hvF/e*R^ for magnetic field (6 — B/Hn is reduced mag- 
netic induction). The free energy density is given in units 
of NqT^ with Nq being the density of states of a normal 
metal. 

To complete the description the selfconsistency equa- 
tions—are required. The gap equation is used it in the 
fornJld 
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The strength of the pairing interaction enters Eq. (Q) 
only implicitly via Tc- We use spherical angles {9, ip) 
to specify v. Then the FS averaging is simply = 
J sin6'd6'd(p/47r. The other selfconsistency equation is for 
supercurrents. It can be disregarded if we concentrate on 
strongly type II superconductors with the GL parameter 
K 3> 1. In this case the internal magnetic field of a su- 
perconductor is constant with high accuracy, of the order 
Removal of this approximation does not invalidate 
the ideas we discuss below, but adds some complexity to 
the formalism and will be considered elsewhere. 

We consider the free energy of |Su superconductor in 
the form introduced by Eilenbergerll^. It is a functional 
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over the order parameter A as well as the quasiclassical 
anomalous Green functions / and / (while g is a depen- 
dent variable by virtue of Eq. (||)). The magnetic term 
(V X A)^ is constant in the approximation we-, use and 
can be omitted. The free energy density readal^l 



F = 



V 



\A\'^lnt + 2ntJ2 
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(1 - g) ^ 2c.„ + ^D/ - Xd*/ 



/A* + /A 
(5) 



where the x integration is performed over the volume V 
of the system. The thermodynamic variables for the free 
energy is temperature and magnetic induction, therefore 
we will discuss t — b plane in what follows. As usual, b is 
related to the external magnetic field by thermodynamic 
arguments. 

Making use of the smallness of A near Hc2{t) line, we 
proceed according to the following strategy. First, we 
solve the Eilenberger equations for / and g perturba- 
tively in A and obtain the gap equation solely in terms 
of A. From that point we deal with this equation only. 
The second step is to work out the Hc2 problem: an 
eigenvalue problem given by the linearized gap equation. 
After Hc2{t) line is known the parameter a/j, which speci- 
fies the distance from this line to an arbitrary (t, b) point, 
can be defined and the expansion for the order parameter 
A itself in terms of au can be constructed. In this paper 
we demonstrate how to obtain the coefficients of the ah 
expansion in general and calculate them, along with the 
free energy, to the lowest order. The Green functions are 
then also known to the same order. 

Step 1. The structure of Eqs. (|l|)-@ suggests the 
expansions: 



9 = 90 + 92 



/ = /i + /3 + 



(6) 



and analogously for /. Subscripts signify the power of 
A to which the term is proportional. The expansions 
can be worked out easily if "inversion" operators P = 
l/(w„ + V • D/2) and P' = - v_- D72) are used 

to solve Eq. (0) for / and Eq. (||) for / . Then, starting 
from the known go = sgn(a;„), we obtain: 



/i=.goPA, fi = goP'A*, 
92 - ~f (pa) (p'A*) , 

(pa) (p'z 
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and so on. The "inverse" operators are conve- 
niently defined using an integral identity: P(v) = 



(see for example dis- 



/„°°dTexp [-t(|c^„|+vD/2J 

cussion in Ref.lli) and P'(v) = P*(-v). Intro- 
ducing creation and annihilation operators by a = 
A {D^ + iDy) /V2, a) = -A {D^ - iDy) /V2 where A eee 
1/a/6 is the magnetic length we arrive at 
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It was assumed that nothing in our problem depends on 
xa, the distance along external magnetic field direction 
(we are dealing with the ground state). 

Step 2. The linearized gap equation is obtained by 
substituting /i from Eq. (0) in Eq. (^. The inte- 
gration produces mi — m2 in the P operator meaning 
that only the number operator a^a and its powers are 
left and, consequently, the eigenfunctions are the Lan- 
dau levels. The highest critical temperature at a given 
external magnetic field is achieved for zero Landau level: 
A(x) ^ ipo{x)n Note that in this case only the first term 
with mi — m2 = in Eq. (p^, which does not depend 
on spatial derivatives, contributes to Hc2- 

Step 3. We are now in a position to introduce the 
actual small parameter ah controlling the problem. For 
this purpose we consider again the term with /i in t 
gap equation Eq. (^ and, following the idea of Ref. 
separate the constant part with m = of the P operator 
from Eq. (10). Then the gap equation can be written as 
follows: 



where 
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and /a, /s, ... have to be taken from Eqs. (0)-(||)j 
etc. As discussed above, the Hc2{t) line is given by 
ah{t,b) — which is identical to the result of Helfand 
and Werthameia. For an arbitrary [t, b) point the pa- 
rameter ah specifies its "closeness" to the Hc2{t) line. It 
is clear that ah can be used to develop an expansion for 
the function A(x) in the vicinity of the Hc2{t) line where 
a/j <C 1 by construction. Eq. ( [ll] ) dictates: 



A = (Ao + a,jAi 



(14) 



with the structure identical to that found by the ah ex- 
pansion of GL theoryO and constitute, most probably, 
a consequence of the fact that the effective electron- 
electron interaction is treated in the BCS theory of su- 
perconductivity in the saddle point approximation. The 
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expansion of Eq. (|TJ) could work well quite far away 
from Hc2(t). To test this assumption we computed the 
function ah{t, b). As seen in the insert in Fig. 1 the nec- 
essary condition for convergence, < 1, is fulfilled in a 
very large portion of the t ~ b plane. 

By construction the operator H defined by Eq. 
does not contain a constant term. This facilitates all 
manipulations with Eq. (pT|). In particular, it is easily 
checked that Aq contains the zero Landau level only since 
TiAo = 0. The equations in the next two orders read: 

HAi = Ao-2TTtY,f [Ao(PAo)(P'A5)J , (15) 

CO „ 

HA2 = Ai+2^tVi [/3{A2Ai}+/5{A5}] . (16) 

The form of the quantities appearing in the integrand of 
Eq. ( p^ is self-explanatory. For example, /3|AqAi} 
is a sum of the expressions given by Eq. (^) with A 
substituted for by Ag twice and by Ai once in all possible 
combination. We do not write down the lengthy exact 
expressions since we do not use them below. We observe 
that the expansion in powers of leads to the mixing 
of orders of the preliminary expansion of Green function 
in terms of A (see Eq. @). Also, the higher order is a 
term in the ah expansion for A, Eq. (p^), the larger is 
the maximal power of nonlinear contributions to it. 

Step 4- We look for solutions for every A^ in the form 
of expansions in the Landau levels' basis {ipm}- 
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In the complete basis of Landau levels there exists an- 
other parameter, besides m, for labeling each function in 
the set. It is analogous to wave vector k. Functions with 
fc 7^ do not contribute to A for equilibrium solutions 
studied in this paper, but of course would be relevant for 
fluctuationscl. The structure of the expansion is such 
that to find the coefficient of ipo in each order the next 
order equation is needed. In particular, the coefficients c 
and cim with m 7^ are found from Eq. (|l^) while the 
coefficients cio and C2m with m are found from Eq. 
(|l6|), and so on. 

Having completed the construction of the ah expan- 
sion we proceed to work it out to the lowest order: to 
calculate the coefficient c (see Eq. (p7|)). We multi- 
ply Eq. (|l^) by on the left, integrate over x and 
use the orthonormality of Landau level basis. Then, the 
term with Ti. vanishes while the first term on the left 
hand side reduces to a constant and we obtain = 
^tT.u=o§FsL^*oPiMH^){P'ro)]- At this stage the 
ideas developed for GL theory in Ref.lla turn out to be 
very effective. In order to lower the free energy the coor- 
dinate dependence of the order parameter A(x) should 
be that of a periodic lattice. This can be parametrized 
by just a few numbers. They are used to minimize the 
free energy calculated first with a lattice of an arbitrary 



shape. Such a minimization procedure is a necessary step 
for the ah expansion too, but we emphasize that this 
nonperturbative step is required only once, in the low- 
est order. We consider Eq. (||), and constraining A and 
/ to satisfy the Euler equation of this functional, Eqs. 
(|l|)-(0), we obtain the equilibrium free energy density: 
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(18) 



We check that, when A is small, Fequu starts from the or- 
der A" (see Eqs. (0)-(|) and note that go = 1 for n > 0). 
The order A^ has vanished, as it should be. Next we in- 
sert the ah expansion for the order parameter A(x) and 

obtain = -c'alfj:Jp^Jjro{P'ro){Pi'o? + 

cc] in the lowest order, which is needed for the minimiza- 
tion. The set {ipm} of lattice shaped Landau levels we 
used is described in detail in Ref.E3. For parameters spec- 
ifying the geometry of the lattice we made a conventional 
choice of p = (02/01) cosa and a = (02/01) sina where 
oi, 02 are edges of the unit cell parallelogram and a is the 
angle in between. In the isotropic case two parameters 
suffice because the global orientation of the lattice is ir- 
relevant while the area of the unit cell is fixed by the flux 
quantization. The further calculations for both the free 
energy Fgquu and the coefficient c contain similar steps. 
After the explicit form of P operator is inserted into Eqs. 
( p^ ) and ( p8[ ) the (p integration and the w„ summation 
can be performed. Handling numerous m summations 
and T integrals allows to present the free energy density, 
in units of NqT^, as 
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(19) 



where the notation Pe = was introduced. Subscript 
"E" stands for Eilenberger. Eq. ( p^ ) can be viewed as a 
generalization of the corresponding result of GL theory, 
which has the same form but with Abrikosov's energy 
parameter Pa = IV'ol^ instead of Pe- Eilenberger's 
energy parameter is given by: 



Pe = 
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Temperature and magnetic induction enter the sums in 
Eq. (|20|) via the combination d = t/y/b only. The coeffi- 



cients Cmim-2 read: 
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The infinite u integral is cut off either by the Gaussian 
or by u/ sinhu. The range of the former depends on d 
because 
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As a result, two different limiting regimes exist for the 
energy parameter: (3E{t,b) ^ if d <C 1 and 

Psitjb) ^ 1/t a d ^ 1. All geometrical information 
about the lattice is contained in matrices /3mim2 — 
'0o'0mi+m2V'mi '01*12 '^hich Can be cast in the form: 



same manner in the entire T ~ H plane. If an approxi- 
mate expression for the Hc2 (T) line is used as a starting 
point, the resulting expansion would have a smaller 
range of applicability. The presented method seems to 
be practical because the actual shape of FS can be con- 
veniently modelled in the framework of the quasiclassical 
approach. 
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where Hm[---] is Hermite polynomial. 

To find the global minimum of FequU from Eq. ( |l9| ) 
we evaluated Pe{p,<^) along the boundary of thcJunda- 
mental domain for these geometrical parametergla. This 
boundary corresponds to all possible two fold symmetric 
lattices and this is the location of extrema. We consid- 
ered (t, b) values belonging the Hc2{t) line and found that 
the minimum of Pe is always at p = 1/2, a = ^3/2, 
that is the hexagonal lattice. Along the Hc2{t) line 
/3£;(l/2, \/3/2) changes only shghtly: from 0.90 at t=l to 
1.23 at t=0 (see Fig.l). Note that it is finite at t = even 
though the P operator, Eq. (|l0|), is poorly defined at this 
temperature. The calculation of /3e along the Hc2{t) line 
produces also the values of Pe for all {t, b) points be- 
low this line because Pe is a function of d = t/\/b only, 
disregarding the 1/t factor common to all terms in Eq. 
(§)• 

The GL results appear from microscopic expressions 
in the limit d oo (or t — > 1 , 6 — > 0) . This limit prPr 
duces only an asymptotic series, not a divergent oneo. 
However, the linear Hc2{t) curve, standard in GL theory 
and obtainable just with the first two terms of the 1/d 
expansion in the equation = (see Eq. (^2|)), is in an 
excellent agreement with the exact Hc2{t) line down to 
as low as t = 0.5. Similar situation turns out to occur 
for the free energy. Our analysis has revealed that even 
at low temperatures the first term with mi = m2 = 
in Eq. (|2^) exceeds the sum over all Landau levels by 
about 15% only (see Fig.l). This interesting fact could 
be one of the reasons why GL theory enjoyed so many 
successes while used very far beyond the domain of its 
strict validity. 

In conclusion, we developed a microscopic descrip- 
tion of type II superconductivity based on the expansion 
in "distance" from the Hc2{T) line, the ah expansion. 
The isotropic case was treated in detail as an example. 
Anisotropic FS and/or pairing interaction would produce 
a different Hc2{T) line. Still, as long as it can be found 
exactly the au expansion is possible to construct in the 
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1-t, GL result for 
Hc2(tj line 
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FIG. 1. Energy parameter /3e for the hexagonal lattice 
along Hc2{t) line presented as a function of temperature: 
dashed line — the contribution of the first term in Eq. ( |2c| ) 
with mi — m2 = 0, solid line — the full sum. Insert: Contour 
plot of the parameter ah as a function of t and b. Magnetic 
induction is normalized to produce db/dt = — 1 at i = 1. 
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